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ABSTRACT 

We present an analysis of deep mid- to far-infrared (MFIR) Spitzer photometric observations 
of the southern 2Jy sample of powerful radio sources (0.05 < z < 0.7), conducting a statistical 
investigation of the links between radio jet, AGN, starburst activity and MFIR properties. This 
is part of an ongoing extensive study of powerful radio galaxies that benefits from both complete 
optical emission line information and a uniquely high detection rate in the far- infrared (far-IR). 
We find tight correlations between the MFIR and [OIII]A5007 emission luminosities, which are 
significantly better than those between MFIR and extended radio luminosities, or between radio 
and [OIII] luminosities. Since [OIII] is a known indicator of intrinsic AGN power, these correla- 
tions confirm AGN illumination of the circum-nuclear dust as the primary heating mechanism for 
the dust producing thermal MFIR emission at both 24 and 70^m. We demonstrate that AGN 
heating is energetically feasible, and identify the narrow line region clouds as the most likely 
location of the cool, far-IR emitting dust. Starbursts make a major contribution to the heating 
of the cool dust in only 15-28% of our targets. 

We also investigate the orientation dependence of the continuum properties, finding that the 
broad- and narrow-line objects in our sample with strong emission lines have similar distribu- 
tions of MFIR luminosities and colours. Therefore our results are entirely consistent with the 
orientation-based unified schemes for powerful radio galaxies. However, the weak line radio galax- 
ies (WLRG) form a separate class of objects with intrinsically low luminosity AGN in which both 
the optical emission lines and the MFIR continuum are weak. 

Subject headings: galaxiesiactive - infrared:galaxies 
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1. Introduction 

In the last few decades space-based infrared ob- 
servatories have provided new opportunities to un- 
derstand the physics of AGN. In particular, mid- 
to far-infrared (MFIR) observations can address 
many outstanding issues relating to the origin of 
the prodigious emission, as well as bringing in- 
sight to questions relating to the unification of ac- 
tive galaxies, the triggering of the activity, and the 
evolution of AGN in general. 

Several efforts have been made to model the 
MFIR spectral energy distributions (SEDs) of 
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radio-loud AGN in order to understand the origin 
of the MFIR emission. The thermal emission from 
warm dust radiating in the mid-infrared (mid-IR) 
is readily explained in model s by AGN heating of 
dust close to the active core dPier fc Krolikl|l992l 
IPier fc Kroliklll993llGranato fc Danesejll994l ). but 
the origin of the thermal far-infrared (far-IR) 
emission is less well understood. This is because 
the far-IR continuum is emitted by cooler dust 
(< 50K), and models that attempt to explain the 
heating of the cool dust in terms of AGN illumina- 
tion of a compact optically thick dust torus have 
struggled to explain how radiation from the AGN 
can heat dust at sufficiently large radii to pro- 
duce emission in the far-IR. Alternative theoreti- 
cal studies that focus on adjusting the geometry 
and size of the torus, or modelling it as a collection 
of discrete dust cloud s surrounding the AGN (e.g. 
Nenkova et alj 20021 van Bemmel fc Dullemondl 



20031 iFritz et al.ll2006l ), can account for the far-IR 
emission in terms of AGN illumination. However, 
the main problem with testing such models is our 
lack of direct knowledge of the spatial distribution 
of the emitting dust. Due to their large distances, 
the circum-nuclear dust structures remain unre- 
solved for the vast majority of AGN. Without 
constraints on the radial distribution of the dust, 
there is a limit to how much the SED modelling 
approach can inform us about the origin of the 
MFIR emission. 

Statistical studies that correlate MFIR, opti- 
cal and radio properties provide a promising al- 
ternative to direct, spatially resolved, studies of 
the dust. This is due to the fact that the MFIR 
continuum is likely to be more isotropic than the 
shorter wavelength UV-optical-near infrared con- 
tinuum, and also because the thermal emission 
from circum-nuclear dust acts as a bolometer for 
the AGN activity. MFIR observations of radio- 
loud AGN are particularly important because the 
extended radio lobe emission from such galaxies is 
generally considered to be emitted isotropically, 
providing an opportunity to select orientation- 
unbiased samples. 

Many previous studies have acknowledged these 
benefits ( Golombek et al.lll988 L Impev fc Gregorini 



19931 IHeckman et alj|l992l Il994 iHes et al.lll 995. 



Haa s et al.l 12004 IShi et alj 120051 IClearv et al 



the infrared observatories, as well as the lack of 
complete, homogeneous samples of powerful ra- 
dio galaxies in the local Universe. It is notable 
that IRAS and ISO detected less than 30% and 
50% respectively of powerful 3C radio galaxies at 
moderate redshifts (z < 0.3). IRAS-based stud- 
ies such as tho se by IHeckman et al. (1994) and 



Hes et al.l (|1995l ) , established putative correlations 



between low frequency radio and MFIR emission 
for radio-loud AGN, indicating a link between the 
AGN power and the MFIR continuum emission. 
However, due to the limited detection rate in the 
MFIR, the heating mechanism for the dust pro- 
ducing the prodigious MFIR emission remained 
uncertain. 

It is now widely accepted that the mid-IR 
(<30/im) continuum is heated by direct AGN il- 
lumination of dust structures close to the AGN 
(e.g. the circum-nuclear torus). However, in light 
of the large scatter of correlations between op- 
tical and far-IR (> 30/im) continuum proper- 
ties, it has been suggested that illumination by 
a starburst component may provide the prin- 
cipal heating mechanis m for the dust produc - 
ing the far-IR emission ( Rowan- Robinson 1995). 
More direct observational evidence for starburst 
heating of the far-IR con t inuum in AGN is pre- 
sented in Schweitzer et al.l ( 20061 ) and Netzer et al 



(l2007t ). )ased on an analysis of starburst-sensitive 
polycyclic aromatic hydrocarbon (PAH) features. 
Finding a correlation between PAH 7.7/mi and 
60/im l uminosity in a sampl e of nearby PG 
quasa rs, (Schweitzer et al. (200(3) and Netzer et al 



2007). However, definitive results were hampered 



in the past by the relatively low sensitivity of 



(|2007l ) attribute this to a link between far-IR lumi- 
nosity and star formation. However, their sample 
is modest, and the PAH star formation signature 
remains undetected in 60% of the objects. There- 
fore the putative correlation between PAH and 
far-IR continuum properties lacks a solid statis- 
tical foundation. While there is no doubt that 
starburst heating of cool dust can account for 
a substantial fraction of the far-IR flux in some 
AGN, the starburst contribution to the far-IR con- 
tinuum in the general population of AGN remains 
uncertain. 

As well as their importance for understanding 
the main dust heating mechanism, MFIR observa- 
tions can also be used to test t he orientation-based 
unified schemes (|Barthell ll989). under the assump- 
tion that the MFIR emission is isotropic. To date, 
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the results from such tests have been ambiguous . 
Early IRAS-bas ed studies ( Heckman et al.l ll994. 



Hes et al.l 119951 ) presented evidence for stronger 



MFIR emission in broad line radio galaxies and 
quasars (BLRG/Q) compared with narrow line ra- 
dio galaxies (NLRG), suggesting that the MFIR is 
not, in fact, isotropic. On the other hand, addi- 



tional studies using ISO data ( Meis enheimer et al 



I2001L lHaas et al]l2004T ) found no evidence for dif- 
ferences between the MFIR luminosities of the 
two optical classes. Unfortunately, these stud- 
ies were hampered by the relatively poor sensi- 
tivity of the IRAS and ISO observatories. Fur- 
ther studies using Spitzer found evidence for a dif- 
ference betw een the MFIR em ission of BLRG/Q 



and NLRG. IShi et all (j2005l ) attribute this dif- 
ference to an isotropic emi s sion at mid-IR wave- 
lenghts, while lCleary et al.l (|20071 ) attribute it to a 
combination of non-thermal contamination of the 
MFIR emission as well as anisotropic emission at 
mid-IR wavelenghts. However, both these results 
were based on samples that were heterogeneous 
and/or incomplete in terms of far-IR detections. 

Many of the past studies of MFIR emission 
from radio galaxies have selected samples based 
on the 3C radio catalogue. Exploitation of this 
catalogue is at present hampered by the lack of 
published high quality optical spectroscopic ob- 
servations for many of the objects, with which 
one can confidently classify and identify possi- 
ble links between MFIR emission, starburst and 
AGN a ctivity. In contrast, the southern 2Jy 
sample (jTadhunter et al.l [1993J) is unique in the 



sense that deep spe ctra have been published for 
the whole sample (TTadhunter et ail Il993l Il998l 
l2002fc iWills et alj[2004 iHolt et alj|2007t See jgj). 
The completeness and availability of deep spectro- 
scopic and radio data make this sample well suited 
to investigating the nature of the MFIR emission 
and testing the unified schemes. Therefore we 
have undertaken a program of deep imaging with 
Spitzcr/MIPS of the 2Jy sample, in order to ad- 
dress the sensitivity problems of previous MFIR 
observat ories leading to inc omplete sample statis- 
tics (see lDicken et al.ll2008l hereafter D08). 

A preliminary analysis of the Spitzer/MIPS 2Ja 



A preliminary analysis 01 the bpitzer/lviLrb zjy 
data set was presented in lTadhunter et aL ( 20071 ) 
and the measured Spitzer MFIR fluxes are pre- 
sented in D08. In this paper we conduct an in- 
depth analysis of the results: sections 2 and 3 



present the sample and the data; section 4 is con- 
cerned with the origin of the MFIR emission; sec- 
tion 5 analyses the far-IR emission and the con- 
tribution of starburst heating to the emission at 
these wavelengths; and section 6 discusses these 
results in the context of the heating mechanism, 
the covering factor of the MFIR-emitting dust, the 
slopes of the correlations and the unified schemes 
for powerful radio galaxies. 

2. The sample 

The sample selected for this study comprises 
a complete sample of all 46 powerful radio galax- 
ies and steep-spectrum quasars (F v oc v~ a ,a\ % > 



0. 5 jj] s elected from the 2 Jy sample of lWall fc Peacock 
(jl985l ) with redshifts 0.05 < z < 0.7, flux den- 
sities 52.7GHz > 2Jy and declinations S < 10°. 
This sample is a comple te, redshift lim ited, sub- 



set of that presented in iTadhunter et alJ (|1993f ) , 



with the addition of PKS 0347+05, which has 
since proved to fulfil the sam e selection criteria 
( di Serego-Alighieri et al.| [T994) . Our selection cri- 
teria cut out all the 16 quasar-like objects in the 
full 2Jy sample that only meet the 2Jy flux crite- 
rion because of the strength of their beamed flat 
spectrum radio core/inner jet components^; the 
remaining objects in the sample are all dominated 
by their extended steep spectrum lobe/hotspot 
emission. Therefore, given that there is no strong 
evidence for beaming and anisotropy in extended 
steep spectrum radio components, our sample is 
unlikely to be significantly biased towards a par- 
ticular orientation of the jets to the line of sight. 
The lower redshift limit has been set to ensure 
that these galaxies are genuinely powerful sources. 
Further discussion of the sample selection can be 



1 In addition to excluding quasars with af'f > 0.5, we also 
excluded the quasars PKS0159-11 and PKS0842-75 on the 
basis that they have relatively strong unresolved radio core 
emission that pushes them above the 2Jy flux limit for the 
sample as a whole. 

2 The only potentially ambiguous cases are PKS0521-36 and 
3C273, which have flat spectra based on their integrated 
radio emission, but extended steep spectrum emission com- 
ponents with flux S2.7GHZ > 2Jy. However, given that the 
extended steep spectrum emission in these objects is con- 
centrated in one-sided jets that are likely to be strongly 
beamed and hence likely to be dominated by non-thermal 
emission, we have decided not to include them in the 
analysis presented in this paper. We find that the inclu- 
sion/exclusion of these objects in the statistical tests makes 
no difference to the main conclusions of the paper. 
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found in D08. 

In the detailed analysis of the spectral energy 
distributions presented in D08, we showed that a 
maximum of 30% of our complete sample have the 
possibility of contamination of their MFIR flux by 
non-thermal synchrotron emission. This is con- 
sistent with sev eral previous studies of samples 
of rad io sources ( Polletta et al. 2000L Clearv et al 



2007), which have indicated that the number of 
objects with possible non-thermal contamination 
of the MFIR is small, and generally confined to 
quasars with flat spectrum radio cores. 
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Table 1 
The Sample 



PKS Other Optical Radio z L 2 4(W/Hz) L 70 (W/Hz) L[om](W) i^ifo ( W / Hz ) Starburst 
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.6 X 
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.0 X 
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WLRG 


FRII 





.116 


1. 


.5 


X 


-i ™23 

10 


1.3 


X 


-i ™23 

10 


2 
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X 


-1 ™33 

10 


4 


.2 x 


10 


No 


0105 — 16 


3C32 


NLRG 


FRII 





.400 


1. 


.8 


X 


10 


< 2.2 


X 


10 


1. 


1 


X 


10 35 


3 


.0 x 


10 


No 


0117-15 


3C38 


NLRG 


FRII 





.565 


3 


.7 


X 


-i ™24 

10 


1.2 


X 


10 25 


7 


.2 


X 


10 35 


7 


.3 x 


10 26 


No 


0213-13 


3C62 


NLRG 


FRII 





.147 


8. 


9 


X 


., „23 

10 


8.2 


X 


* ^23 

10 


5 
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X 


10 


4 


.4 x 


- ~25 
10 


No 


0235-19 
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BLRG 


FRII 





.620 


5. 


.5 


X 


-i ™24 

10 


7.1 


X 


. „24 

10 


8 


.(> 


X 


10 35 


9 


.8 x 


10 26 


No 


0252 — 71 




NLRG 


CSS 





.566 


1. 


.7 


X 


-i ™24 

10 


< 5.5 


X 


. „24 

10 


6 


.4 


X 


-. ™34 

10 


9 


.7 x 


10 26 


No 


0347+05 




BLRG 


FRII 





.339 


8. 





X 


-i ™23 

10 


7.1 


X 


i ™24 

10 





.1 


X 


1 ™33 

10 


2 


.0 x 


10 26 


U 


0349-27 




NLRG 


FRII 





.066 


4 


.2 


X 


10 22 


2.0 


X 


10 23 


5 


.3 


X 


10 33 


9 


.2 x 


. ..24 

10 


U 


0404+03 


3C105 


NLRG 


FRII 





.089 


2 


.6 


X 


10 23 


6.0 


X 


10 23 


1. 


.3 


X 


10 34 


2 


.0 x 


10 25 


No 


0409-75 




NLRG 


FRII 





.693 


2 


.3 


X 


., „24 

10 


1.7 


X 


10 25 


5 


.(> 


X 


., „34 

10 


3 


.6 x 


- ~27 

10 


SB 


0442—28 




NLRG 


FRII 





.147 


5. 


1 


X 


-i ™23 

10 


7.2 


X 


i ™23 

10 


3 


.1 


X 


i ™34 

10 


5 


.6 x 


10 25 


No 


0620-52 




WLRG 


FRI 





.051 


1. 


3 


X 


., „22 

10 


1.4 


X 


., „23 

10 


< 1 


.1 


X 


- ^32 

10 


3 


.3 x 


.. ~24 

10 


SB 


0625-35 


OH-342 


WLRG 


FRI 





.055 


7. 


.6 


X 


-i ™22 

10 


1.4 


X 


i ™23 

10 


1 


.3 


X 


1 ™33 

10 


(i 


.5 x 


-i ^24 
10 


No 


0625-53 




WLRG 


FRII 





.054 


5. 


.5 


X 


10 


< 3.4 


X 


10 22 


< 4 


.9 


X 


10 32 


4 


.7 x 


. ..24 

10 


No 


0806-10 


3C195 


NLRG 


FRII 





.110 


3 


.4 


X 


i ^24 

10 


0.4 


X 


i ™24 

10 


2 


.(> 


X 


10 35 


2 


.2 x 


10 25 


No 


0859-25 




NLRG 


FRII 





.305 


!) 


.2 


X 


10 23 


8.3 


X 


10 23 


4 


.2 


X 


i ™34 

10 


2 


.4 x 


10 26 


No 


0915-11 


Hydra A 


WLRG 


FRI 





.054 


2 


.9 


X 


i ™22 
10 


3.8 


X 


. „23 

10 


1 


.3 


X 


i ™33 

10 


4 


.2 x 


10 25 


SB 


0945+07 


3C227 


BLRG 


FRII 





.086 


3 


.3 


X 


„23 

10 


1.3 


X 


., „23 

10 


3 


.5 


X 


- ^34 

10 


2 


.1 X 


- ~25 

10 


U 


1136-13 




Q 


FRII 





.554 


6 


.1 


X 


i ™24 
10 


1.1 


X 


10 25 


2 


.4 


X 


10 36 


8 


.9 x 


10 26 


U 


1151-34 




Q 


CSS 





.258 


1 


.5 


x 


10 24 


4.7 


X 


10 24 


1 


.2 


x 


10 35 


2 


.3 x 


10 26 


U 


1306-09 




NLRG 


CSS 





.464 


2 


.0 


X 


10 24 


9.1 


X 


10 24 


6 


.3 


X 


10 34 


5 


.9 x 


10 26 


U 


1355-41 




Q 


FRII 





.313 


6 


.0 


X 


10 24 


7.5 


X 


10 24 


3 


.4 


X 


10 35 


2 


.0 x 


10 26 


u 


1547-79 




BLRG 


FRII 





.483 


2 


.9 


X 


10 24 


7.0 


X 


10 24 


1 


.2 


X 


10 36 


5 


.1 X 


10 26 


u 


1559+02 


3C327 


NLRG 


FRII 





.104 


2 


.8 


X 


10 24 


5.5 


X 


10 24 


7 


.8 


X 


10 34 


3 


.4 x 


10 25 


No 


1602+01 


3C327.1 


BLRG 


FRII 





.462 


2 


2 


X 


10 24 


3.5 


X 


10 24 


2 


.9 


X 


10 35 


3 


.4 x 


10 26 


No 


1648+05 


Here A 


WLRG 


FRI 





.154 


6. 


.6 


X 


10 22 


< 6.2 


X 


10 23 


2 


.0 


X 


10 33 


3 


.7 x 


10 26 


No 


1733-56 




BLRG 


FRII 





.098 


3 


.3 


X 


10 23 


1.7 


X 


10 24 


2 


.9 


X 


10 34 


3 


.5 x 


10 25 


U 


1814-63 




NLRG 


CSS 





.063 


2 


.5 


X 


10 23 


5.9 


X 


10 23 


1. 


9 


X 


10 33 


1 


.4 x 


10 25 


U 


1839-48 




WLRG 


FRI 





.112 


4. 


.5 


X 


10 22 


1.6 


X 


10 23 


< 1 


.0 


X 


10 32 


1 


.8 x 


10 25 


No 


1932-46 




BLRG 


FRII 





.231 


2 


1 


X 


10 23 


1.5 


X 


10 24 


1. 


1 


X 


10 35 


2 


.4 x 


10 26 


SB 


1934-63 




NLRG 


CSS 





.183 


6. 


2 


X 


10 23 


7.1 


X 


10 23 


5. 


2 


X 


10 34 


2 


.6 X 


10 26 


No 


1938-15 




BLRG 


FRII 





.452 


2. 


.3 


X 


10 24 


6.6 


X 


10 24 


3 


.4 


X 


10 35 


9 


.2 x 


10 26 


No 


1949+02 


3C403 


NLRG 


FRII 





.059 


6 


.9 


X 


10 23 


1.2 


X 


10 24 


3 


2 


X 


10 34 


8 


.6 x 


10 24 


No 


1954-55 




WLRG 


FRI 





.060 


1. 





X 


10 22 


3.3 


X 


10 22 


< 4 


.5 


X 


10 31 


6 


.1 X 


10 24 


No 


2135-14 




Q 


FRII 





.200 


4 


.5 


X 


10 24 


4.9 


X 


10 24 


5 


.9 


X 


10 35 


6 


.5 x 


10 25 


U 


2135-20 


OX-258 


BLRG 


CSS 





.635 


6 


2 


X 


10 24 


5.3 


X 


10 25 


6 


.4 


X 


10 35 


1 


.1 X 


10 27 


SB 


2211-17 


3C444 


WLRG 


FRII 





.153 


1. 


.8 


X 


10 22 


< 3.4 


X 


10 23 


1. 


1 


X 


10 33 


6 


.2 x 


10 25 


No 


2221-02 


3C445 


BLRG 


FRII 





.057 


7. 


.4 


X 


10 23 


5.9 


X 


10 23 


7. 


.6 


X 


10 34 


7 


.6 x 


10 24 


No 


2250-41 




NLRG 


FRII 





.310 


1. 


.4 


X 


10 24 


2.7 


X 


10 24 


2 


.2 


X 


10 35 


1 


.7 x 


10 26 


No 


2314+03 


3C459 


NLRG 


FRII 





.220 


4. 


.0 


X 


10 24 


4.0 


X 


10 25 


7 


.2 


X 


10 34 


8 


.1 X 


10 25 


SB 


2356-61 




NLRG 


FRII 





.096 


4. 


.0 


X 


10 23 


7.3 


X 


10 23 


4 


.0 


X 


10 34 


4 


.7 x 


10 25 


No 



Note. — Tabled Column 3 definitions; Q - quasar, BLRG - broad line radio galaxy, NLRG - narrow line radio galaxy, 
WLRG - weak line radio galaxy. Column 4 definitions; FRI & FRII - FanarofT- Riley class 1 and 2 respectively, CSS - com- 
pact s teep spectrum, C/J - co re/jet. Column 8 [OIII]A5007 luminosities calculated from flux presented in Tadhuntcr ct al. 
Jl993h and lWills eTaLl £2004). Column 10 gives the 15-17Ghz monochromatic core luminosities from the fluxes presented 
in DOS. Column 11 gives information about whether a spectroscopic young stellar population (YSP) has been detected at 
optical wavelengths indicating the possibility of starburst activity. SB - YSP detected, No - No YSP, U - uncertain YSP 
component (references for YSP/starburst: |Tadhunter~et"a l. 2002 Wills ct al. 2004, 2008; Holt et al.ll2007l) . Positions for 
the objects can be found in D08. 



5 



3. The data 

We have made deep Spitzer/MIPS observations 
(24, 70 and 160/im) of our sample, as well as 
complementary high frequency radio observations 
(15 to 22 GHz) with the ATCA and the VLA. 
Full details of the observations and reduction can 
be found in D08, along with the MFIR, radio 
fluxes and spectral energy distributions for the en- 
tire sample. We detect 100% of our sample at 
24/im, 90% at 70/xm and 33% in the lower sen- 
sitivity 160/xm band. This is by far the best de- 
tection rate for MFIR observations of a sample of 
an intermediate-redshift AGN published to date. 
Additionally, the high frequency radio data, along 
with data from the literature, have enabled us to 
detect the radio cores in «70% of our complete 
sample. We utilize these data to calculate the 
orientation-sensitive R parameter (see D08 and 



m 
line 



In Table [T] we present the 24 and 70/im 
monochromatic luminosities for the complete sam- 
ple, along with 5 GHz radio and [OIII] A5007 emis- 
sion line luminosities^. The analysis presented 
here differs from the prelim inary results presented 
adhunter et al.l (|2007 ) in that the emission 
luminosities have now been corrected for 
Galactic extinction, using E(B-V) reddening val- 
ues obtained from the NASA/IPAC Extragalactic 
Database (NED), along with the param eterised 
Galactic extinction law of lHowarthl (|l983l) . In Ta- 
ble [T] we also present the most up-to-date optical 
classification for each of the sources, with objects 
classified as narrow-line radio galaxies (NLRG), 
broad- line radio galaxies or quasars (BLRG/Q), 
and weak-line radio galaxies (WLRGjS 

In Figure [1] we plot 5 GHz radio luminosity 
against both redshift (top) and radio spectral in- 
dex (bottom) for the complete sample. The effect 
of the flux limit of the sample is clearly visible in 



3 In order to calculate luminosities we used H = 
71kms- 1 Mpc- 1 ,fl m = 0.27 and S7 A = 0.73 along with 
spectral indices derived from the F(70)/F(24) flux ratios 

4 WLRGs are sometimes known as low-excitation galaxies 
but we prefer to label them as WLRG since the excitation 
(or, more accurately, ionization) of the emission line gas 
is not necessarily related to the AGN luminosity. Indeed, 
there exist examples of AGN with relatively high [OIII] lu- 
minosities, but emission line ratios reflecting a low ioniza- 
tion state. WLRG are defined as having EW([OIII])<10A 
( Tad hunter et al.lll998t) . 
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Fig. 1. — Distribution of 5 GHz vs. redshift to- 
tal radio luminosity (top); high frequency radio 
spectral index (o^f) vs - 5 GHz radio luminosity 
(bottom) for the 2Jy sample. 

the upper plot in the form of the tight correlation 
between monochromatic radio luminosity and red- 
shift. It is noteworthy that our complete sample 
covers three orders of magnitude in radio luminos- 
ity. The plot of the high frequency radio spectral 
index a|;f against radio luminosity allows us to in- 
vestigate any bias in our steep spectrum selection 
method. Based on the similar distributions of the 
different optical classes in a^'j over the full 3 or- 
ders of magnitude covered by our sample, we are 
confident that the steep spectrum selection has left 
us with no strong bias towards BLRG /Q or NLRG 
classification. 

4. MFIR emission and the central engine 



Ma ny previous stu d ies (e.g 



19941 iHes et al.l Il995l IClearv et all 



Heckman et al 
20071) 



have 

found evidence that MFIR emission power is cor- 
related with measured total radio power, albeit 
based on highly incomplete samples. This result 
is in line with a model in which the MFIR emis- 
sion and radio jet power are both strongly linked 
through the physics of the central engine. How- 
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Fig. 2.— Luminosity correlation plots: (a) L 2 4/im vs. £[orn]A5007; ( b ) vs. £[oni]A5007; (c) ^24^™ vs. L 5 gh z ; (d) L 70lim vs. 

isGHz- The cross in the bottom right corner of plots (a) and (b) represents the maximum uncertainties in [OIII] and MFIR luminosity 
measurements, demonstrating that the scatter in the plots is real and not a consequence of observational uncertainties. Red circles indicate 
BLRG/Q objects, green squares NLRG and blue triangles WLRG. 



ever, radio emission depends on factors in addition 
to the intrinsic power of the AGN. For example, 
the properties of the local ISM interacting with the 
radio plasma are expected to have a major impact 
on the conversion of jet mechanical power int o ra- 
dio luminosity (e.g. Barthel fc Arnaud 1996). In 
addition, the low frequency radio emission may be 
emitted by extended structures that are far from 
the regions emitting the MFIR continuum. There- 
fore, a change in the intrinsic power of the AGN 
may take a significant amount of time to be re- 
flected in the emission of the extended radio lobes. 

An alternative to comparisons based on low fre- 
quency radio emission is to use optical emission 
line luminosity. The AGN-photoionized narrow- 
line region (NLR) is emitted on a smaller scale 
(<5kpc) than the extended radio lobe emission in 
most radio galaxies. Therefore the [OIII]A5007 
emission line is likely to provide a good indica- 
tion of the intrinsic power of the i l luminating AGN 



e.g. j Rawlings & Sa undersill991l iTadhunter et al 



19981 ISimpsonlll998l) . 



4.1. Comparisons with MFIR luminosity 

We now utilize our complete sample to address 
the question of how the MFIR continuum depends 
on the AGN power, with the benefits of not only 
completeness in the MFIR flux measurements, but 
also accurate [OIII]A5007 emission line luminosi- 
ties for the entire sample. 

In Figures^ (a) and (b) we plot 24//m and 70/im 
monochromatic luminosities (£24/™, Lro^m) against 
[OIII] luminosities (£[om])> while in Figures [2] (c) 
and (d) we plot the 24/mi and 70/im luminosities 
against 5GHz monochromatic radio luminosities 
(^5GHz)- It is clear from Figures 2 (a) and 2(b) 
that the MFIR luminosities are strongly correlated 
with Lrorm. As discussed in the preceding sub- 
section, the [OIIIJA5007 luminosity is expected to 
provide a good indication of intrinsic AGN power. 
Therefore, we propose that the 24/zm and 70/im 
emission are also intrinsically linked to the power 
of the active core and, consequently, that the dust 
producing the MFIR emission is likely to be heated 
directly by the central AGN. 

It is also notable from Figure [2]J a) that L24/im 
shows a much tighter correlation with the L[om] 
than it does with the Lsghz- This result is ex- 
pected, because we believe that both the 24/im 



luminosity and the [OIII] luminosity are strongly 
linked to the power of the active core via AGN 
illumination of the emission line clouds and dust 
structures close to the nucleus. The radio power, 
however, is dependent on additional factors such 
as the nature of the ISM in the halo of the host 
galaxy into which the jets and lobes expand. 

Considering next the differences between the 
24/im and 70/im correlations, it is apparent that 
there is more scatter in the L-jQ^ m vs Lroni] corre- 
lation. LjQ^ m also remains better correlated with 
L[oui] than Lsghz, but the difference is not as ev- 
ident as at 24/tm. Our interpretation of the cause 
of this additional scatter at 70/im lies in the origin 
of the far-IR emission, as we will discuss in detail 
in Sj5] However, it is noteworthy that this scatter 
cannot be due to the reduced sensitivity of Spitzcr 
at longer far-IR wavelengths because our observa- 
tional errors are small compared with the overall 
scatter in the distribution (see error bars in plots 
(a) and (b)). Therefore, there must be a physical 
cause for the increased scatter in this plot. 

4.2. Rank correlation statistics 

Prior to interpreting our results in full, it is im- 
portant to investigate the significance of the corre- 
lations in Figure 2. We have therefore calculated 
the Spearman rank correlation coefficient for the 
four correlations. The test was undertaken using 
a redshift limited sub-sample with z>0.06 in or- 
der to remove the majority of WLRGs with upper 
limits in [OIII]. In addition the remaining object 
with an upper limit in [OIII] (PKS1839-48) was 
also excluded^. In order to quantify the effects 
that the four remaining upper limits in 70/im have 
on this test, we chose to handle the correlations 
that included 70/im data in the following way: we 
created replacements for each of the four 70/im 
upper limits by randomly selecting a 70/im/24/tm 
flux ratio value from the distribution of measured 
70/im/24/im flux ratio values for the sample, and 
then multiplying this by the measured 24/im flux 
to create a new 70/im flux estimate. The 70/im 
estimates were then converted to luminosities and 
included in the rank correlation test. This pro- 



5 The exclusion or inclusion of the [OIII] upper limit value for 
this object has been investigated, and we find that it makes 
no significant difference to the statistical results presented 
in Table [2] 
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Table 2 
Statistical Analysis 



Rank Correlation r s 



(1) L 24 vs L[oni] 


0.88 


(2) L 70 vs L[oni] 


0.76 


(3) L 2 4 vs L 5GHz 


0.54 


(4) L 70 vs L 5GH , 


0.63 


(5) L5GH2 vs z 


0.93 


(6) Lsghz vs L[oiii] 


0.57 


Partial Rank Correlation with z 


(7) L 2 4 vs L[oiii] 


0.77 


(8) L 70 vs L[oiii] 


0.50 



Note.— Tabic [2] Result of var- 
ious Spearman rank correlation 
statistics. Values of < r s < 1 
are given for each test , where a 
value close to 1 is highly significant. 
This test was undertaken with a z 
limited sample z>0.06, to remove 
most of the objects with upper lim- 
its in [OIII]. In addition the object 
PKS1839-48 which has upper limits 
in [OIII] was also removed, leaving 
38 objects from our complete sam- 
ple. Sec P~2l for discussion of the 
70/^m upper limits. 



cess was repeated 1000 times, and the median of 
the correlation coefficient for those cycles is pre- 
sented in Table [2] for correlations involving 70/im 
luminosities (i.e. rows 2, 4 and 8). 

The results show a high level of significance for 
all the correlations using a two tailed test, in the 
sense that we can reject the null hypothesis that 
the variables are uncorrelated at a >99.9% level of 
significance. The rank correlation statistics, pre- 
sented in Tabled clearly show that the Lronii vs. 
^24/jm and Lroiii] vs. £70^™ correlations are more 
significant than those between the radio luminos- 
ity and MFIR and [OIII] luminosities. 

For a flux-limited sample such as that con- 
sidered here, it is natural to consider the possi- 
bility that the correlations between MFIR, radio 
and [OIII] luminosities might not be intrinsic, but 
rather arise through their mutual dependence on 
redshift. For example, given the strong correlation 
between radio power and redshift induced by the 
radio flux limit (see Figure 1), and also the corre- 
lation between emission line luminosity and radio 
power, it is possible that a correlation between the 
L[Oin] an d the MFIR emission luminosity might 
arise because the MFIR luminosities are indepen- 



dently correlated with z (e.g. because of genuine 
redshift evolution rather than being intrinsically 
correlated with £[om])- The second part of Table 
[5] (rows 7 and 8), shows the results of a partial rank 
correlation test carried out in order to investigate 
whether the correlations could be a result of a de- 
pendence on the third variable z. In both cases 
we still find that the null hypothesis that the vari- 
ables are unrelated can be rejected at a >99.5% 
level of significance. This demonstrates that both 
^24/jm and I/70/im are intrinsically correlated with 
£[Oin]- 

4.3. LIRG and ULIRG comparison 

In order to place radio galaxies in the context of 
other sources of prodigious MFIR emission in the 
local universe, it is interesting to compare them 
with the luminous and ultra-luminous infrared 



galax ies (LIRGs and ULRIGs: ISanders fc Mirabel 
119961 ). Figure [3] again shows the L^o^n vs -^[0111] 
correlation for the entire sample, but in this case 
we also indicate the typical ranges that LIRGs 
and ULIRGs would occupy on the diagram. To 
calculate these ranges we follow the definitions of 
Sanders fc MirabelT ([19961 ): that LIRGs have in- 
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Fig. 3. — Plot of LiOfim vs L\oiii] showing the 
comparison with LIRG and ULIRG luminosities. 
Approximate ranges for LIRGs and ULIRGs are 
indicated by the horizontal lines, where the dashed 
red lines, represent the luminosity limits of typi- 
cal LIRGs and ULIRGS for a warm SED model, 
and the dot dashed blue lines represent the lim- 
its of LIRGs and ULIRGs for the cool SED model 
(see H4.3p . Stars represent objects with optical ev- 
idence for recent star formation activity (see jj37 



tegrated far-IR luminosities 10 n Lo > Lfir > 
1O 12 L , while ULIRGs have 1O 12 L > L FIR > 
10 13 L©. The ranges of 70/im monochromatic lu- 
minosities corresponding to these integrated lumi- 
nosities have been estimated based on two assump- 
tions about the MFIR SEDs of the sources. In 
each case the lower line has been calculated by 
assuming the SED of the radio galaxy 3C327 - 
taken to be representative of a source with rela- 
tively warm MFIR colours, whereas the upper line 
has been calculated by assuming the SED of the 
ULIRG Arp220 — representative of a source with 
cool MFIR colours. Considering these limits we 
find that 28-48% of our sample would be classified 
as LIRGs, and <7% as ULIRGS. 

4.4. Testing the unified schemes 

The MFIR luminosities can also be used to test 
the orientation-based unified schemes for powerful 
radio galaxies (regardless of the emission mecha- 
nism). In particular, if the distributions of the 
observed MFIR luminosities are similar for differ- 
ent optical classes of objects (NLRG, BLRG/Q, 



WLRG) this is consistent with, but does not prove, 
that they are part of the same parent population. 

Evidence from inv e stigations using IRAS (e.g. 
Heckman et al.l Il994l iHes et all Il995l ) suggested 
that the MFIR luminosities of narrow line ra- 
dio galaxies (NLRG) are lower by up to an or- 
der of magnitude in comparison to broad line 
radio galaxies and quasars (BLRG/Q These 
results have been supp orted by some studies 
van Bemmel et all 2000h and rejected by o thers 



(jMeisenheimer et al]|200ll lHaas et al. 2004), but 



all the previous studies were based on incomplete 
samples. Therefore, it remains uncertain whether 
genuine differences exist between the MFIR prop- 
erties of broad and narrow line radio-loud AGN. 

Any difference between the two types of radio 
galaxies at shorter mid-IR wavelengths (<30/im), 
might be explained in terms of extinction by a 
dusty torus, provided that there is significant op- 
tical depth in the torus at such wavelengths. How- 
ever, it would remain challenging to explain dif- 
ferences in longer wavelength emission in terms of 
such obscuration. Alternatively, a difference be- 
tween the two classes of objects could be due to 
non-thermal beamed emission contaminating the 
thermal MFIR, which is expected to be stronger 
in the BLRG/Q. 

Using our complete sample we can investigate 
whether there are significant differences between 
the MFIR properties of broad- and narrow-line ob- 
jects in the 2Jy sample. 

A visual inspection of Figures [2] (a) and (b) 
reveals no evidence that the BLRG / Q have higher 
luminosities than the NLRG at 24 or 70/xm, and 
we have shown that no more than 24% of the sam- 
ple have a possibility of contaminati on from non- 



therm al beamed components (see iDicken et al 
(|2008h ). Therefore we conclude that the MFIR 



emission in our sample is most likely to be emit- 
ted isotropically, at least down to an observed 
wavelength of 24/im (rest wavelength 14-23/im, 
depending on redshift). 



3 In the following, unless otherwise specified, we consider 
BLRG and quasars as a single class. This is justified on the 
basis that, among objects with broad permitted emission 
lines, there is a continuous range in luminosities between 
objects classed as BLRG and quasars. 
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Fig. 4. — Plot of the orientation-sensitive R vs. 24/im luminosity (left) and 70/im luminosity (right), where 
R is defined as S cor e/{Stot - S cor e)- 



4.5. The R parameter 

We can further examine whether the MFIR con- 
tinuum is isotropic by investigating any links be- 
tween MFIR luminosities and the relative bright- 
ness of the radio core. If there is any anisotropy 
in the MFIR continuum, we would expect to see 
enhanced emission in those objects with strong 
cores, which have axes orientated closer to the 
line of sight. Using the new high frequency ra- 
dio observations presented in D08, as well as data 
from the literature, we can investigate the de- 
pendence of MFIR properties on the orientation- 
sensitive, core dominance parameter R, defined as 
Score/ {Stot — S core ). Tables of R values estimated 
at 5GHz for the complete sample are presented in 
D08. 

In Figured] we plot the 24 and 70/im luminosi- 
ties against R. It is clear from these plots that, on 
average, the B LRG / Q have the high est R values 



in our sample (jMorganti et alj Il997f ) , consistent 



with unified schemes which predict that BLRG/Q 
have jet axes that are pointing closer to the line of 
sight, leading to stronger beamed core/jet emis- 
sion. However, not all BLRG/Q are highly core 
dominated. In fact, the population of such objects 
shows a wide range of R parameter values. 

In Figure [5] we plot R against the 70/im/24/im 
infrared flux ratio. Plotting similar parameters, 



Shi et al.l (|2005l ) found an anti-correlation in the 
sense that warmer IR colour corresponds to a 
smaller R parameter. In addition, they found 
some evidence that different optical classes lie in 
different regions of the plot, with the BLRG/Q 
tending to the highest R, and warmest colours. 

From a visual inspection of Figure [5] there is 
little evid ence for a strong anti-correlation as pre- 
sented in I Shi et al.l (|2005l) . The range of infrared 
colours suggests that the MFIR emission is rela- 
tively isotropic. However, it is worth noting that 
many of the objects in the lower region of the dia- 
gram, forming the group at low 70/im/24/xm and 
high R, are BLRG/Q. 

4.6. Weak line radio galaxies 

It is not only the distributions of NLRG and 
BLRG/Q that are of interest in the correlation 
plots presented in Figured The figure also pro- 
vides clues to the nature of WLRG. Although 
many of the WLRG in our sample have been classi- 
fied as FRI radio sources, a significant subset hav e 
FRII morphologies ( see Tadhunter et al.l Il998h . 
Previous studies (e.g. ICao fc Rawling sl l2004h have 
suggested that obscuration could provide an ex- 
planation for the differences between WLRG and 
radio galaxies with strong emission lines. 

In Figures [2] (a) and (b) it is striking that 
the WLRG lie at the lowest luminosity end in 
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Fig. 5. — Plot of the orientation sensitive R pa- 
rameter against MFIR colour, represented by the 
70/24/zm flux ratio. 

both emission line and MFIR continuum luminos- 
ity. If the MFIR emission is primarily caused 
by AGN heating of the dust, then this demon- 
strates that WLRG contain intrinsically weak 
AGN. Therefore, WLRG cannot be accommo- 
dated in the simplest unified schemes for radio 
galaxies, which posit that the differences be- 
tween all classes of radio-loud AGN are solely 
due to anisotropy and orientation. It is evi- 
dent that, although obscuration may contribute 
at some level to the weakness of the emission lines 
in these objects, the idea that WLRG contain 
powerful quasar-like AGN, heavily obscured by 
dust, is certainly not consistent with our data. 
Hardcastle et all (120071 ) suggest that the differ- 
ences between WLRG and radio galaxies with 
strong emission lines (NLRG/BLRG/Q) may be 
accounted for by different accretion modes of the 
AGN, with the WLRG accreting hot ISM at a rel- 
atively low Eddington ratio, and NLRG/BLRG/Q 
accreting cold ISM at a higher Eddington ratio. 

Figures (c) and (d) clearly show that, despite 
being weak in both emission line and MFIR lu- 
minosity, the WLRG overlap with the NLRG and 
BLRG/Q in terms of radio luminosity. This sug- 
gest that other factors play an important role in 
boosting the radio powers of WLRG. In particular, 
as noted above, the radio luminosity is not neces- 
sarily a good measure of intrinsic AGN may power, 
and the properties of the local ISM have a large 
impact on the radio luminosity for a given intrinsic 



jet power. In this case we would expect that, for a 
given MFIR luminosity, the WLRG with powerful 
radio emission should be found in relatively dense 
cluster enviro nments. Some evidence for this effect 
is presented in Barthel fc Arnaud ( 19961 ). however, 
their result is based on a limited sample. Future 
deep optical imaging observations of the 2Jy sam- 
ple will provide a direct indication of the signifi- 
cance of such environmental effects. 

5. Origin of the Far-Infrared Emission 

As discussed in the nature of the heating 
mechanism for the cool dust radiating at 70/im 
(the far-IR) is not well understood. However, it 
is widely accepted that the cool dust is likely to 
be heated by starbursts and/or by AGN illumi- 
nation. Based on the tight correlation between 
£[om] an d Lmfir shown in Figure [2j we have 
concluded that the most probable heating mech- 
anism for the dust emitting at both mid- and 
far-IR wavelengths is direct illumination by the 
AGN (Q. However, it is of course naive to rule 
out starburst heating as a contributor to far-IR 
emission altogether, specifically because morpho- 
logical evidence suggests that at least some pow- 
erful radio galaxies are triggered in maj or, gas 



rich galaxy mergers (|Heckman et al.lll986l ). Such 



mergers are predict ed to be associated with power- 
ful starbursts (e.g. Idi Matteo etaD (|2005h ). The 
starburst- AGN connection is also important for in- 
terpreting sub-millimetre observations in the con- 
text of the star formation history of radio-loud 
A GN at high redsh i ft. Som e studies, such as that 
bv lArchibald et all (|200lh . assume that the cold 
dust responsible for the sub- mm emission is heated 
entirely by starbursts. Therefore to comprehen- 
sively understand the MFIR emission from pow- 
erful radio galaxies it is vital to investigate the 
contribution of starburst heating to the observed 
far-IR fluxes, and possible links between star for- 
mation and AGN activity. 

5.1. Evidence for starburst heating of the 
far-IR continuum 

Careful spectral synthesis modellin g of the high 
quality optical spectra for our sample (ITadhunter et al 
I2002L IWills et al.ll2004 [20081 iHolt et alJl2007h has ~ 
allowed us to identify the objects that show clear 
evidence for recent star formation activity at op- 
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Fig. 6. — Plots showing the variation of Lmfir with £[om] f° r ^ ne complete sample at 24 and 70/im, with 
those objects identified as having young stellar populations marked with separate symbols (stars). Also 
displayed is a supplementary sample of all known radio galaxies with spectroscopically-confirmed evidence 
for star formation activity at optical wavelengths, the properties of which are detailed in Table El The 
regression line has been fitted only to the data from the original sample with z > 0.06, in order to avoid 
most of the objects with upper limits in their [OIII], and also excluding the 7 starburst objects. 



tical wavelengths (see Table [T]). Therefore, in 
Figure [6] we are able to plot the L24^ m and I^o/im 
data against £[oin] once again, but this time high- 
lighting those seven objects that we have estab- 
lished to have unambiguous evidence for starburst 
activity at optical wavelengths (PKS0023— 26, 
PKS0409-75, PKS0620-52, PKS0915-11 (3C218), 
PKS1932-46, PKS2135-20, PKS2314+03 (3C459)). 
As we have seen in Sj4j there are good correla- 
tions between both L^m, Lj ^ m an( i ^[Oiii]> 
but there is significantly more scatter in the 
70/xm correlation , consi stent with the findings of 
iRowan-Robins^r] (l995). As first described m the 



prelim inary results presented in iTadhunter et al 
(|2007h . we believe this scatter is in part due to 



the boosting by starburst heating of the cool far- 
IR emitting dust that is not as significant for the 
warm mid-IR 24/im dust emission. 

To demonstrate the degree of starburst boost- 
ing above the main correlations we have plotted re- 
gression lines on both plots in Figure [5] The lines 
shown are the bisector of a linear least squares fit 
of x on y and y on x. As for the sample used in the 
Spearman rank statistics, in calculating this fit we 
include only those objects with z > 0.06, in order 



to avoid most of the objects with upper limits in 
their [OIII] emission. In addition, the 7 starburst 
objects and one more object (PKS1839— 48) with 
upper an limit in [OIII] were also removed, leaving 
a total sample of 33 objects for the fit. 

From a visual inspection of Figure Elb) it is 
clear that the majority of those objects we identify 
as having evidence for optical star formation activ- 
ity are displaced, showing enhanced far-IR lumi- 
nosities relative to the regression line in the 70/im 
plot; this enhancement is not apparent for most 
of the optical starburst objects in the 24/im plot. 
Therefore, it appears that there must be some 
mechanism boosting the far-IR flux in starburst 
compared to non-starburst objects. Analysing the 
displacements of starburst objects we find that 
their far-IR fluxes are boosted by a factor of four 
on average relative to the remainder of the sample, 
with a maximum boosting factor of 20 in the case 
of 3C459. 

In addition to our complete sample, we have 
also plotted a supplementary sample of starburst 
radio galaxies, taken from the literature on Fig- 
ure El These supplementary sources represent all 
known radio-loud AGN that show spectroscopic 
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evidence for star formation activity at optical 
wavelengths, apart from the similar objects in the 
2Jy sample. The properties of these sources are 
presented in Table [31 and they are plotted as open 
stars in Figure O It is clear that these supple- 
mentary objects also lie well above the regression 
line fitted to our data in the 70/im plot. A visual 
inspection of the supplementary starburst sample 
in Figure [6] also shows a tendency for increased 24 
and 70/im luminosities, relative to the majority of 
the sample, at low L[om] emission. This is not 
surprising given that, at low [OIII] luminosities, 
the AGN itself is intrinsically weak and a modest 
amount of ongoing star formation can boost the 
MFIR fluxes well above the regression line. 

Due to the high level of completeness of our 
sample in terms of MFIR detections, our data can 
also be used to test the statistical significance of 
any differences between the starburst and non- 
starburst populations in Figure [5] 

First, for our complete sample we have consid- 
ered a one dimensional Kolmogorov- Smirnoff (K- 
S) two sample test, comparing the vertical dis- 
placements from our fitted regr ession line in the 
L jQfm vs. Lroiii] plot (see also iTadhunter et al 



2007). In this case we find that we can reject 
the null hypothesis that the starburst and non- 
starburst samples are drawn from the same par- 
ent population at a better than 1% level of signif- 
icance. We have made the same test, this time in- 
cluding both the 2 Jy starburst and the supplemen- 
tary starburst objects, and investigating whether 
they are significantly displaced relative to the re- 
gression line for the non-starburst 2Jy sample. In 
that case we find that we can reject the null hy- 
pothesis that the starburst and non-starburst sam- 
ples are drawn from the same parent population 
at a level of significance of better than 0.1%. 

Secondly, we consider the significance of the 
displacement between the starburst and non- 
starburst samples using a two dimensional K 



S test (Peacock 19831 Fasano fc Franceschini 



19871 ). The method we apply here is the gen- 
eralisation of the 2D K-S test developed by 
Fasano fc Franceschinil (1987). Again for the 



Lrn^m vs i[oni] plot, initially considering the 
starburst objects against the non-starburst ob- 
jects (in our complete 2Jy sample), we find that 
we can reject the null hypothesis that the star- 
burst and non-starburst samples are drawn from 



the same parent population at a level of signifi- 
cance of better than 5%. We find exactly the same 
significance level when including the supplemen- 
tary starburst objects. In addition, applying the 
same statistical test to the £[om] vs - plot 
we find no differences between the starburst and 
non-starburst samples, as expected given that the 
24/im luminosity is less likely to be effected by 
starburst heating. 

There are three notable outliers from the cor- 
relation, which are labelled on Figure [6{b). These 
sources add to the increased scatter in the far-IR, 
and we describe them in more detail below. 

• PKS0347+05. This object lies behind the 
Galactic plane and, despite the relatively 
large correction of the [OIII] flux for Galac- 
tic extinction, it still lies above our fitted 
regression line. In a recent analysis of near- 
IR images for the source it was found that 
there are, in fact, two AGN within our MFIR 
photometric aperture for this object. If 
both AGN radiate at MFIR wavelengths this 
could perhaps explain some of the apparent 
excess relative to the correlation. Addition- 
ally, since PKS0347+05 is a BLRG, its AGN 
is likely to outshine any starburst signatures 
at optical wavelengths. Therefore, it is pos- 
sible that this object could also have star 
formation activity that has not so far been 
detected at optical wavelengths. 

• PKS0945+07. This object lies below the 
correlation, and is an example of an object 
with an inverted MFIR spectrum that ap- 
parently lacks a c o ol dust component (e.g. 



2001) 



Mile v et all 11984 Ivan Bemmel fc Bar thel 



PKS1306-09. This ob ject was discussed in 
Tadhunter et al.l (|2002h . and shows marginal 



polarization at optical wavelengths, which 
could be due to scattered light or a non- 
thermal optical continuum component. In 
addition, the extrapolation of the radio com- 
ponent in the SED (see D08) is consistent 
with non-thermal contamination of the ther- 
mal MFIR continuum. So far, the optical 
spectra are inconclusive regarding the pres- 
ence of a young stellar component in this 
object, but it is difficult to entirely rule out 
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the presence of such a component. 



Table 3 

The Supplementary Starburst Sample 



Object z L 24 (W/Hz) L 70 (W/Hz) £ [0 ni](W) L 5GH z(W/Hz) SB rof [OIII] ref 



3C48 


0.367 


3 


.0 


X 


10 25 


1 


.6 


X 


10 26 


3 


.0 > 


< 10 36 


1 


.6 


X 


10 26 


6 


4 


3C213.1 


0.194 


5 


.1 


X 


10 24 a 


5 


.0 


X 


10 24 a 


5 


.6 > 


< 10 33 


3 


.1 


X 


10 25 


18 


11 


3C236 


0.101 


2 


.0 


X 


10 23 


5 


.7 


X 


10 23 


1 


.8 > 


< 10 33 


1 


.7 


X 


10 25 


1,5 


2 


3C285 


0.079 


9 


3 


X 


10 23 a 


1 


.7 


X 


10 24 a 


3 


.1 > 


< 10 33 a 


4. 


.1 


X 


10 24 


1,2 


2 


3C293 


0.045 


9 


.2 


X 


10 22 


6 


.7 


X 


10 23 


1 


.3 > 


; 10 32 


4. 





X 


10 24 


7 


1 


3C305 


0.042 


9 


.4 


X 


10 22 a 


5 


.1 


X 


10 23 a 


1 


.4 > 


( 10 34 a 


1 


.7 


X 


10 24 


7 


3 


3C321 


0.096 


3 


.3 


X 


10 24 


1 


.0 


X 


10 25 


8 


.8 > 


( 10 34 


1 


2 


X 


10 25 


1,3 


2 


3C433 


0.102 


2 


.1 


X 


10 24 


3 


.5 


X 


10 24 


2 


.2 > 


( 10 34 


4 


2 


X 


10 25 


1.4 


3 


B2 0648+27 


0.041 


1 


.5 


X 


10 24 a 


4 


.(> 


X 


10 24 a 


1 


.2 > 


( 10 34 a 


8. 


.6 


X 


10 22 


12 


7 


B2 0722+30 


0.019 


1 


.6 


X 


10 23 a 


1 


.1 


X 


10 24 a 


5 


.6 > 


; 10 31 a 


2 


4 


X 


10 22 


13,14 


8 


Con A 


3.4Mpc 


8 


.5 


X 


10 21 a 


1 


.2 


X 


10 22 a 


5 


.0 > 


; 10 33 


2 


4 


X 


10 23 


15,16 


9,10 


PKS0131-36 


0.030 


1 


.5 


X 


10 23 a 


1 


.6 


X 


10 24 a 


2 


.5 > 


; 10 32 a 


3 


.5 


X 


10 24 


1 


5 


PKS0320-37 


0.006 


7 


.8 


X 


10 21 


1 


.5 


X 


10 23 


3 


.0 > 


< 10 31 


2 


.5 


X 


10 24 


8,9 


5 


PKS0453-20 


0.035 


1 


.2 


X 


10 23 a 


8 


.8 


X 


10 23 a 


6 


.9 > 


< 10 31 a 


2 


.1 


X 


10 24 


11 


5 


PKS1345+12 


0.122 


1 


.2 


X 


10 25 a 


3 


.4 


X 


10 25 a 


9 


.0 > 


( 10 34 a 


4. 


.6 


X 


10 25 


1,10 


6 


PKS1549-79 


0.152 


1 


.2 


X 


10 25 a 


2 


.7 


X 


10 25 a 


8 


.9 > 


( 10 34 a 


1 


.1 


X 


10 26 


17 
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Note. — Table [3] supplementary sample comprising of all radio-loud objects from the literature with spectro- 
scopic evidence for recent star formation activity at optical wavelengths. All infrared fluxes were obtained 
from NED apart from the object 3C305 for which Spitzcr/MIPS fluxes were extracted from the archive im- 
ages. 24 a nd 70pm lumino s ities f or 3C 213.1 are derived from IRAS upper limits. SB refer ences: ( 1 ) [Holt_e^_alj 
12007'! . (2) lAretxaga et al.l 1120011') . (31 ITadhunter et al.l 119961) . (4) IWills et al.l 12002). (51 IQ'Dea et all l|200lF 
f6) ICanalizo fe Stockton! 1200l|)~(7 1 ITadhunter et all 12005ft .78) IGoudfrooii et all 120011) , (9 )IKuntschnerl iboool) , 
(lOl lRodriguez Zaurin et al.l 12007ft . (lllwills et a.1,1 <l200-tfl . (1 2iEmonts et al.l l|20061) . ( 13~l JEmont A I l|2006ft . (141 
Emonts et al. in preperation, (151 IPeng et al. (2002 1. (161 IPeng et al.1 ^2004 | ). (171 [Holt et al l 120071) , (18) 

IWills et al] 120081) . [OIII] flux references: (1) lEmonts et al.l 120051), (21 | Saunders et all 119891). (3) R obin- 
son PhD thesis, 2001 Unive rstiy of Sheffield, (41 IChatzichristou et all 11993 1, (51 ITadhunter et al.l (jl993D. (6 1 

IGelderman & Whittle] 1199411 . (7) | Emonts et al l 12006ft . (81 Emonts private communication. (91 ISimPBonl jl998fl . 
f 10) IStorchi-Bergmann et al.l 119971) , f 11) IHoltl 120051) . 

^derived from IRAS measurements at 24/jm and 60pm. 
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5.2. Correlation slope statistics 

In the calculation of our fitted regression line 
plotted in Figure [51 we have chosen not to include 
those objects identified as having evidence for a 
star burst component. Using a bootstrap tech- 
nique we have investigated the uncertainty in the 
slopes of our fit to these correlations and tested the 
effect that the removal of the 70/im upper limits 
and starburst objects would have on the regression 
lines. 

In order to investigate the uncertainties, a sam- 
ple of N data points (x and y points from Figure[6|) 
were labelled and then drawn at random to create 
a group of N replacements for the original sam- 
ple. We then recalculated the slope and repeated 
this process 1000 times. Provided the data points 
are independent, the distribution of the slopes esti- 
mated in the bootstrap trials provide an indication 
of the uncertainty in the estimated slopes. 

The data fitted for this bootstrap are identi- 
cal to those used in the Spearman rank correla- 
tion statistics ( 34. 2j) i.e., with redshifts limited to 
z > 0.06 in order to avoid most of the objects 
with upper limits in their [OIII] luminosity and 
PKS1839— 48. Here we investigate the effects of 
including starburst objects and the upper limits 
in the sample drawn for the N data replacements, 
in order to quantify their effect on the fitted slopes. 

The starburst objects are simply included or 
excluded when calculating the slope however, ex- 
actly as the Spearman rank test (<0T2j), the 70/im 
upper limits were included by randomly choosing 
a ratio of 70/im/24/im to derive new 70/im upper 
limits. This process was repeated for each of the 
1000 bootstrap trials. The results of our analysis 
are presented in Table [4] 

Column 2 of Table [4] contains the values of the 
slopes we have fitted to the data in Figure [6] (the 
bisector of a linear least squares fit of x on y 
and y on x) as well as the slope values calculated 
with starburst objects excluded. Columns 3 and 4 
present the mean value for 1000 slopes of the sam- 
ple, created using a bootstrap technique, column 
3 includes the 70/im upper limits and column 4 is 
for the sample bootstrapped without upper limits. 

Considering, in the first instance, the results 
for the L24/im vs. £[om] correlation, we find that 
the addition/removal of the starburst objects and 
the objects with upper limits in 70/im has little 



effect on the resulting mean slope values, where 
all values are consistent within the estimated lex 
uncertainties. 

Secondly, looking at the results in Table Q] for 
the LyQ^ m vs. i[oin] slopes, it evident, that al- 
though the upper limits do not significantly affect 
the mean regression line slopes and associated un- 
certainties, the introduction of the starburst sam- 
ple does. This further supports our hypothesis 
that the starburst objects are responsible for much 
of the scatter in the Lroiii] vs. LjQ^ m correla- 
tion. Excluding the optical starburst objects, the 
slopes and associated uncertainties of the £70^™ 
vs £[0111] correlation are entirely consistent with 
those of the hi^m, vs £[0111] correlation, reinforc- 
ing the idea that the dust radiating at both mid- 
IR and far-IR wavelengths has a common heating 
mechanism: AGN illumination. 

Given that the bootstrap sample in column 3 
is identical to that used to plot the regression line 
in Figure [6l these statistics map the uncertainty 
in our slope fit. For the sample with upper limits 
but not including the starburst objects, we calcu- 
late that the uncertainties in our regression line 
slopes are 0.05 and 0.07 for the 24/im and 70/im 
correlations respectively, from the standard devi- 
ations of the bootstrap slopes. 

6. Discussion 

6.1. Heating mechanism 

On the basis of our observational results, we 
have concluded that the MFIR continuum emit- 
ted by powerful radio galaxies is predominantly a 
consequence of AGN heating of the circum-nuclear 
dust. However, far-IR emission from dust heated 
by starbursts is apparent in the minority of the 
objects that fall above the main correlations in- 
volving 70/im. 

The link between far-IR bright radio sources 
and optical evidence for starb urst activity has 
been noted in previous studies. iTadhunter et al" 
(|2002j) found that the two objects with the 
strongest evidence for optical starburst activity 
in their sample are als o the most luminous i n 
the far-IR. Subseque ntlv lWills et al. I <|2002Ll2004h . 
while investigating the UV excess in radio galax- 
ies, also found that those objects with the best 
evidence for young stars have significantly larger 
far-IR luminosities than the rest of their sample. 
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Table 4 
Bootstrap Analysis 



Correlation plot 


Bisector 


with 70^111 UL 


without 70/^m UL 


24/xm vs [OIII] with starburst 
24/im vs [OIII] without starburst 


0.78 
0.75 


0.79 ± 0.05 
0.75 ± 0.05 


0.76 ± 0.08 
0.71 ± 0.07 




70/im vs [OIII] with starburst 
70/im vs [OIII] without starburst 


0.83 
0.71 


0.88 ± 0.11 
0.72 ± 0.07 


0.93 ± 0.13 
0.78 ± 0.09 



Note. — Tabled Slopes for the complete sample excluding objects with z > 0.06 and 
PKS1839 — 48. The uncertainties on the values presented in columns 3 and 4 are the 
standard deviations of the slopes derived from the 1000 bootstrap cycles. 



The Spitzer results for this 2Jy sample strongly 
reinforce these earlier results with a much larger 
sample. 
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Fig. 7. — Histograms of MFIR colours 
(70/im/24/xm flux ratio) for the non-starburst 
sample, our starburst sample and the supplemen- 
tary starburst sample. The black line shows the 
dividing line between what we define as cool and 
warm colours (70/im/24/im = 5). 



Cont rary to our results, previous studies such as 

those by Rowan- Robinson dl995 ). van Bemmel fc Barthell 



(|2006h and iNetzer et ahl (|2007h analyzed PAH 
starburst tracers and MFIR emission. Although 
una ble to detect PAH featu res in 60% their sam- 
ple, Schweitzer et all ( 2006I ) suggest that star for- 
mation heating most likely accounts for at least 
30% of the far-IR emission of their sample of PG 
quasars, finding a correlation between PAH lumi- 
nosity at 7.7^tm and 60/im luminosity similar to 
that of starburst dominated ULIRGs. Although 



Schweitzer et al.l ([2006) provide some of the best 
empirical evidence for a starburst/ AGN connec- 
tion, albeit with many upper limits both in the 
far-IR and PAH emission, it is not clear that PAH 
features are exclusively starb urst signatures. In 



deed, studies such as that bv lSiebenmorgen et al 



(|200ll ); lvan Bemmel fc Dullemondl (|2003l ). put for 
ward starburst heating as the primary heating 
mechanism of the dust producing far-IR emission 
in powerful AGN, albeit without direct evidence. 
This has encouraged some investigators to search 
for more direct observational ev idence of a star- 
burst/AGN link. For example, ISchweitzer et al 



(|2004l ) have suggested that it may be possible to 
produce PAH emission from AGN heating alone if 
the dust is situated at sufficiently large radii from 
the AGN. 

We have shown that the far-IR continuum lumi- 
nosities (I^o^m) of the 2Jy radio sources are corre- 
lated with their [OIII] luminosities [L^oni]) with 
a similar slope to the correlation between h^^m, 
and Lyoni]- The similarity of the two correlations 
points towards a common mechanism for produc- 
ing both the mid and far-IR emission. If we as- 
sume, as is generally accepted, that the continuum 
emission at 24/im is heated almost exclusively by 
the AGN then, in order for the cool dust emitting 
in the far-IR to be heated predominantly by star- 
bursts rather than the AGN, a remarkable degree 
of coordination between AGN and starburst activ- 
ity would be implied. Although we cannot entirely 
rule out this possibility, we regard it as unlikely. 

Empirically, starburst galaxies are known to be 
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associated with cool MFIR colours, as expected 
in the case of illumination of extended dust struc- 
tures by spatially distributed star forming regions. 
In Figure [7] we present histograms showing the 
70/im/24/im MFIR colours of the non-starburst, 
starburst and supplementary samples separately. 
It is immediately apparent from the plot that 
all the 2Jy starburst objects have cool colours 
(70/i77i/24^m > 5), whereas all but two of the 
non-starburst sample have much warmer colours 
(70/im/24^m < 5). The non-starburst object in 
Figure [7] with the coolest colour is PKS 0347+05; 
as discussed above, it is possible that this object 
may contain as yet undetected star formation at 
optical wavelengths. Moreover, one of the non- 
starburst objects whose starburst status is un- 
certain - PKS1306-09 - also has relatively cool 
colours (70^m/24/im = 4.7). 

The data presented in Figure[6]can also be used 
to consider links between AGN and starburst ac- 
tivity. As discussed in ij5.ll at low emission line 
luminosities (£roiii] < 10 33 W) it is evident that 
L\oui] is only weakly correlated with the far-IR 
luminosity. However, it is notable that all the 
objects with large 70^m luminosities (L^Q^ m > 
10 25 W/Hz) — many of which we have identified 
with starburst heating of the cool dust — also 
have large [OIII] luminosities (£[om] > 10 35 W). 
This suggest a weak link between starburst and 
AGN activity in the sense that only the most pow- 
erful radio-loud AGN are associated with power- 
ful, ULIRG-like starbursts. However, the reverse 
is not true since many of the highest emission 
line luminosity objects — some of which are lu- 
minous enough at MFIR wavelengths to be clas- 
sified as LIRGs — do not show evidence for star- 
burst activity, suggesting that powerful AGN are 
not always accompanied by massive star forma- 
tion activity. Similar trends have been found 
in a study of nearby type 2 S eyfert galaxies by 
Gonzalez Delgado et al. ( 200ll ). who show that, 
for a given [OIII] luminosity, there are two pop- 
ulations of Seyferts: one with optical evidence for 
prodigious recent star formation, large far-IR lu- 
minosities and cool MFIR colours; the other with 
little evidence for enhanced star formation activ- 
ity. However, in the case of the Seyferts, a higher 
proportion (>50%) of objects show evidence for 
recent star formation activity compared with our 
sample of radio galaxies. 



We can also estimate the rate of energetically 
significant starburst activity in our sample by 
considering the main optical and infrared indica- 
tors of starbursts: 7 (15%) of the objects in the 
2Jy sample show unambiguous spectroscopic ev- 
idence for recent star formation activity at opti- 
cal wavelengths; 9 (20%) have cool MFIR colours 
(L 7 o f _ im /L 2 4^n > 5); 12 (26%) of the objects lie 
more than 0.3 dex («factorx2) about the regres- 
sion line in the Lfo^m vs I/[om] correlation in Fig- 
ure 6; and 13 objects (28%) show at least one of 
these indicators. Therefore we estimate that the 
proportion of powerful radio-loud AGN showing 
evidence for energetically significant recent star 
formation activity is in the range 15-28%. 

As first described in iTadhunter et al.l (|2007l ). 



this brings us to a key result: given the lack of 
a correlation between starburst and AGN activ- 
ity, and the fact that only a minority of objects 
in our sample show any evidence for recent star 
formation activity, it is unlikely that all powerful 
radio galaxies are triggered at the peaks of ma- 
jor gas rich mergers. This conclusion is strength- 
ened by the fact that not all powerful radio galax- 
ies show tidal tails or oth er evidence for recent 
galaxy merger activity (e.g Hcckma n et al. I I1986L 



Tadhunter et al.lll989l . iMcLure et al.lll999h 



6.2. Energetic feasibility 

Having found preliminary evidence, based on 
the correlations between MFIR and [OIII] lumi- 
nosities, that the MFIR emitting dust is predom- 
inantly heated by AGN, it is important to assess 
whether this heating mechanism is energetically 
feasible. We start by assuming a simple model 
in which the far-IR continuum, mid-IR continuum 
and [OIII] emission lines are produced by AGN 
illumination of structures with covering factors 
Cfi r , C m ir and C n i r respectively. We make no 
assumptions about the radial distribution of dust, 
although the far-IR emitting dust must be situ- 
ated at larger radial distances from the AGN than 
the mid-IR emitting dust in order to produce its 
cooler temperature. 

Based on simple recombination theory, the to- 
tal H[3 luminosity (Lup), generated in the narrow 
line region (NLR) of an AGN, is related to the 
ionising luminosity (Li on ) by 
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->H/3- 



(hu) 



a fl ,, 

ton ^eff 



hu 



c: 



(7 + l)L bo lCr, 



_lil.5xl0 14 Hz : 
llxl0 13 //2 



where (hu)i on is the mean ionising photon energy, 
huH/3 is the energy of an H(i photon, affj is the 



is 



B 

From 



total case B recombination coefheient and a 
the effective H0 recomb i nation coefficient 
the results of lElvis et all (|l994h the bolometric lu- 
minosity Lboi is related to the ionising luminosity 
by L bo i 



3.1Li on , therefore 



L H p = 0.32L 



bol 



hv H p 



l H/3 



{hv)i 



At high luminosities and ionization parame- 
ters, L[om] ~ VLLjip. In addition, fits to the 
line ratios of nearby radio galaxies are consis- 
tent with (hu) ion — 6.2 x 10 -U J, corresponding 
to an ionisin g continuum shape w ith = 1.5 
(F v oc v~P) (iRobinson et al.lll987). A lso, given 
that agfl/a e // = 0.09 (|Osterbrocklll989t ). we find: 



and we can write a similar expression for the Kfi r . 
Given our assumption of a simple power-law shape 
for the MF1R SED, the MFIR spectral index is 
directly calculated from the MFIR colour: 

_ ^0.9(^70/^24) 
7 ~ /o<?(24/70) 

For our complete sample we find a median 
MFIR colour of L 7nfim / L 2 A tim = 2.1, leading to 
Jmedian ~ 0.7. Therefore the MFIR monochro- 
matic luminosities can be related to the AGN 
bolometric luminosities by: 

Lmir — 3 X 10 LholCmir^ 



and 



Lfi r — 1.25 X 10 ^LboiCfi r V 



0.7 



L 



[OIII] 



= 2.2 x 10- 2 L bol C n 



lr- 



(1) 



This result can be related to the MFIR luminosi- 
ties by assuming that the structures producing the 
MFIR absorb the short wavelength radiation and 
re-radiate it in the MFIR, then L^° l ir = C m i r L bo i, 
and similarly L b jf r — Cf ir L bo i. Defining the wave- 
length range of the mid-IR as 2 - 30^m, and that 
of the far-IR as 30 - 100/im, the mid-IR and far-IR 
bolometric luminosities are defined as: 



1.5xl0 14 ffz 



Lmir — / L v d,V 
'lxl0 13 ffz 



and 



t bol 



lxl0 1J Hz 



L,,du. 



3xl0 12 Hz 



Representing L v as a power law L v = Kv 1 , we 
can then write 



L 



bol 



1.5xl0 14 Hz 



Ku^dv 



lxl0 13 flz 



and 



By substituting equation (1) into these expres- 
sions we find: 



= 1 X 1U Lii 



and 



L(7Qnm) = 9 x 10 



'[OIII] 



'[OIII] 



Cfir 



(2) 
(3) 



It is likely that the major uncertainty in C m i r /C n i r 
is associated with the assumed L bo i/Li on ra- 
tio that may be uncertain by a factor of ~ 2 
(lElvis et al.lll994fh 



In Figure[5]we have plotted the L? Qlim vs £[01111 
and £24/im vs £[oiirj correlations with predictions 
from our calculations for a range of C m j r /C„j r and 
Cfi r IC n i r covering factor ratios. This allows us to 
investigate whether the MFIR emission can be ex- 
plained solely in terms of AGN illumination of the 
dust. Note that our calculations do not take into 
account the fact that the Lmfir vs £[oiii] rela- 
tionships are non- linear (discussed further in H6.3j) : 
the calculated ratio slopes are most relevant to the 
high luminosity ends of the correlations. Also, in 
the calculation so far we have assumed that the 
NLR, the mid-IR and far-IR emitting regions rep- 
resent discrete structures at different locations in 
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io 35 io 36 io 37 io io~ io~ 

L [ONI] 5007 (W) L [Olll] 5007 (W) 



Fig. 8. — Plots of the £[oiii] VS.L24/J (left) and L[oni] vs.Lyo^ (right) for the complete 2Jy sample with lines 
showing the range of possible ratios of covering factors required to explain the trends in terms of the AGN 
illumination model. C m i r /C n i r in plot (a) and Cfi r /C n i r in plot (b). Note our calculation of the covering 
factors does not seek to explain the non-linear slope (discussed further in <j6.3[) . Also the ranges of covering 
factor ratios plotted are different in plots (a) and (b). The arrow in the upper left hand corners of the plots 
show the effect that correcting for an intrinsic extinction of E(B-V) = 0.3 would have on the [OIII] emission 
line luminosities. 



each galaxy, whereas, in fact, it is more likely that 
there exists a continuous distribution from one re- 
gion to the next, or an overlap between the regions. 

Firstly, Figure [S^a) shows that C mir /C n i r « 5- 
25 would account for the mid-IR luminosities of 
the majority of the objects in our sample. How- 
ever, a value of C rn i r /C n i r « 12.5 provides the 
best description of the high luminosity end of the 
main correlation (from Figure [6]). Therefore NLR 
clouds cannot be responsible for the bulk of the 
mid-IR emission of the radio galaxies in our sam- 
ple, a result consistent with detailed models for the 



mid-IR emission from NLR clouds ( Groves et al 



2006) . We further note that much of the scatter in 
the £24^771 vs ^[oiii] correlation can be explained 
in terms of variations in C mir /C n i r . 

The absolute covering factor of the dusty torus 
structure emitting the mid-IR continuum can 
be estimated by assuming that the NLR has a 
covering factor in the range C n w ~ 0.0 2-0.08 
(jNetzer fc Laorlll993l rMaiolino et all l200ll ). For 

Cmir/C n l r RJ 12.5 this gives 0.25 < C m ir < 1- 

A large C m i T is entirely feasible in the context of 
obscuration by the putative dusty torus required 
by the orientation-based unified schemes (| Bart he! 



1989), since a torus with a typical opening half 
angle in the range 45° — 60° would have a covering 
factor of 0.5-0.7. This calculation also indicates 
that the covering factor of the NLR must be small 
(Cnir < 0.08) if we require C m ir <1. Therefore, it 
is evident that the 24/im continuum emission from 
the warm dust component can be feasibly powered 
by AGN heating alone for most of the objects; the 
relative covering factor of the mid-IR emitting 
structures is consistent with that expected for a 
circum-nuclear torus, which is not itself expected 
to radiate significant forbidden line radiation due 
to high densities and temperatures. 

Secondly, it is clear from Figure [8jb) that the 
ratio of covering factors required for the emitting 
dust structures is much less in the far-IR than at 
mid-IR wavelengths: typically Cf i r /C n i r ~ 1 — 2 
for the main non-starburst sample, correspond- 
ing to Cfi r = 0.02 — 0.16. Such covering fac- 
tors are feasible if cool dust is associated with 
the outer parts of the torus, a kpc-scale dust lane 
suc h as those found in many nearby radio galax- 
ies (|de Koff et al.ll2000h . or indeed the NLR itself. 
In the latter context we note that, if we correct 
the emission line luminosities for the reddening 
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typical of the NLR of radio ga l axies (0.3<E(B- 



V)<1 0: ITadhunter et all 11994 iRobinson et &L 
200(1 lRobinsorJl200lH 7 l this will shift many of the 
points in Figure [SJb) close to the Cfi r /C n ir = 1 
line. Moreover, the NLR in many nearby ra- 
dio galaxies and PG quasars are extended on 
sufficiently large spatial scales (~lkpc or larger: 
Bennert et aDbool IPrivon et all l2008h that the 
dust associated with them would be heated by the 
AGN to the relatively cool temperatures required 
to produce the far-IR radiation. Therefore, rather 
than being separate entities, it is entirely plausi- 
ble that the NLR clouds are the dusty structures 
emitting the far-IR continuum. 

We have again marked the starburst objects 
from our sample in Figure [8] Many of these star- 
burst objects require Cfi r /C n i r w 10 — 50. If 
we attempt to explain the far-IR emission in such 
objects in terms of AGN illumination alone, this 
raises the question of how the AGN could illumi- 
nate the far-IR emitting regions without produc- 
ing substantial emission line radiation from the gas 
associated with the dusty structures (i.e. lead- 
ing to Cfi r /C n i r = 1). One explanation is that 
the C n ir is underestimated, and Cfi r /C n i r over- 
estimated, due to dust extinction of the NLR. 
However, this would require more substantial red- 
dening (E(B-V) > 1) than is supported by most 
observations of the NLR. Therefore the objects 
that lie above the Cfi r /C n i r = 1 line are likely to 
have their far-IR flux boosted by starburst heat- 
ing or non-thermal emission. Incidentally, this in- 
cludes PKS0347+05 and PKS1306-09, which we 
discussed in §5. li as possible starburst candidates. 

6.3. Correlation slopes 

The fitted regression lines from Table |4] (0.75 
±0.05 at and 0.72 ±0.07 at 70^m) are 

plotted on Figure [6] Taking into account the 
bootstrap uncertainties, the slopes of the cor- 
relations a re significantly non-lin ear. It is no- 
table that iMaiolino et al.l (|2007l ) find a similar 
slope (0.82±0.02) between continuum luminosity 
at 5 100 A and infrared luminosity at 6.7/xm for 
a sample of 25 high luminosity QSOs. We now 
consider possible explanations for the non-linear 



7 These reddening estimates are based on Balmer line ratio 
measurements obtained after careful subtraction of the un- 
derlying stellar continuum. 



slopes. 

First we consider explanations related to the 
emission line physics of the NLR clouds. The 
most basic AGN photoionization models assume 
that the emission-line regions are photoionised by 
the AGN and that the properties of the emission 
line region, such as density and covering factor, do 
not change significantly with AGN power. Con- 
sideri ng such a model f or a s imple optically thick 



slab, ITadhunter et al.l (|1998I ) found that at high 



luminosities, the relation between L 



[OIII] 



and 



ising luminosity can be represented as a power- 
law L ion oc (i[ iii]) Q , where a = 0.8 ± 0.1. As- 
suming that the properties of the torus do not 
change with luminosity (but see below) and that 
the MFIR continuum luminosity is directly pro- 
portional to L{, i and Li on , we therefore expect 
r0.8±o.i This result is consis- 



Lmfir oc £[0111] 



oc L 



[OIII] 



tent with the slopes of our Lmfj r vs £[dttt] cor- 



relatio ns as well as that found by IMaiolino et al 
(J2J307J). However, based on several detailed stud 



ies o f the emission line spectra of radio gala xies 
(e.g. IRobinson et~aL1 |2000L rTavlor et al.ll2003r ). it 
is likely that the single slab photoionization mod- 
els are over-simplistic, and that multi-component 
models for the NLR are more appropriate. In this 
case the relationship between L[oni] and Li on is 
less clear. Moreover, investigations of the rela- 
tionship between £[oni] a- n d the continuum lumi- 
nosities in large samples of AGN provide direct 
observational evidence for steeper slopes in the 
relationship between Lhni and £[0771], albeit with 
larger scatter (a > 1 iNetzer et afll2004 Eiblh . 

Alternatively, we can also consider whether the 
non-linear slopes are consistent with the basic re- 
ceding torus model that has been suggested to ex- 
plain the variation in the BLRG±Q/NLRG frac- 
tion with luminosity in the context of the unified 
schemes (lLawrencel(l991 , Hill et al1ll996l[simpson 
19981 l2005h . This solution is attractive because 



it is widely accepted that, if the mid-IR emitting 
dust lies close to the AGN, the inner radius of 
the torus at which the dust is sublimated due to 
illumination by the AGN must depend on AGN 
power. From simple thermal equilibrium argu- 
ments we can relate the distance of the illuminat- 
ing AGN from the inner face of the torus (r) to the 
bolometric luminosity (Lf,oi) and the sublimation 
temperature of the dust (T su b) as follows: 
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r oc 



^bol 
sub 



1/2 



(4) 



The area of the inner face of the obscuring torus 
can be written as Anrh where h is the thickness of 
the torus from its mid-plane to its top surface. 
The covering fraction of the torus is then: 



C t , 



Airrh 
Aitr 2 



h 

r 



(5) 



assuming all the radiation absorbed by the face 
of the torus is re-emitted by dust, and h is fixed, 
then the relationship between emission in the mid- 
IR at 24/im emitted by dust close to the AGN and 
the bolometric luminosity can be written as: 



P24fim OC — Lbol oc L 



1/2 
bol ■ 



having substituted equation (5) into equation (4). 
Finally, assuming that Lbol oc (Lromi) we ob- 
tain: 



Plium OC (£[OIII] 



a/2 



In the context of the receding torus model we 
now consider three cases with different relation- 
ships between I/roin] and LboV- 



In, I 



1. The relationship between £[om] cind L. 
is linear with a=l, and therefore the non- 
linear slopes of i[oni] vsLmfir. correlations 
are entirely due to the receding torus. In 
this case P^Ajim oc (iroiii]) which is too 
shallow to account for the measured corre- 
lation slopes, even when considering the un- 
certainties on the measured slope we have 
estimat ed from the bootst rap analysis (± 
0.05). iNetzer et aD (|2004T ) found evidence 
for such a linear relationship between Lpm] 
and Lagn, from the lack of variation of the 
equivalent width of [OIII] as a function of 
continuum luminosity. However, there are 
many upper limits for the high luminosity 
objects in their sample. 

2. The correlation between L[om] and Lbol is 
not linear and has a slope a < 1, as pre- 
dict ed by simple AGN phot oionization mod- 
els i Tadhunter et al. 199& ). If a=0.8, as 



discussed above, then P2Aam oc (iroini) 4 - 
This value is also significantly shallower than 
measured from our data. 

3. The correlation between £[oiii] a nd Lagn 
is non-linear and a > 1. In t his co ntext 
it is notable that Netzer et al. ( 20061 ) find 
evidence for a « 1.7 ± 0.07, based on an 
analysis of X-ray and emission line prop- 
erties of luminous AGN. This then gives 

„ / n0.85±0.04 . . 

P24 fJ ,m oc [L [Qui] J , which is consis- 

tent with the measured slopes within the un- 
certainties. 

Option 1, considering only a receding torus 
model, cannot satisfactorily explain our results. 
ISimpsonl(|2005l) also shows that, assuming isotropic 
[OIII] emission, the fraction of BLRG/Q increases 
with AGN luminosity, but with a shallower slope 
than predicted by a pure receding torus model. 
He proposes an increase in torus height with lumi- 
nosity to account for the discrepancy. Although 
this may not be th e only ex p lanati on for the ob- 
served trends, the ISimpsonl (|2005l ) result shows 
that there remain questions about the application 
of the most basic receding torus model. 

In this study we have shown that the slope of 
the LiQftm vs i[oni] correlation is similar to that of 
the Li^m vs i[oin] correlation. If the far-IR radi- 
ation is emitted by extended cool dust structures 
it is clear that these must lie far from the AGN for 
the typical luminosities of our sources (> O.lkpc). 
In this case it is unlikely that the dust structures 
will recede in the same way as predicted for the 
torus. Therefore, while a receding torus might pro- 
vide a plausible explanation for the non-linearity 
of the L24um vs i[oni] correlation (for a > 1) such 
an explanation is unlikely to apply to the L-jQ^ m 
vs £[oni] correlation. 

6.4. Unification 

Given the completeness of our sample, the simi- 
larities between the MFIR properties of the NLRG 
and BLRG/Q classes provide excellent support for 
the orientation-based unified schemes for radio- 
loud AGN. In addition, our results provide an im- 
portant constraint for understanding the structure 
and optical depth of the obscuring torus. The 
question then remains: how do we reconcile our 
results with previous studies that did find evidence 
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for significant obscuration at MFfR wavelengths? 
Apart from sample selection and incompleteness 
(see JQ), there are several factors that may have 
biased previous studies and led to the finding of 
a difference between the two classes of radio-loud 
AGN: 

• Classification of BLRG objects: The 

deep spectra available for all our sample ob- 
jects allow us to make secure classifications 
of objects into NLRG, BLRG/Q and WLRG 
categories. Past studies may have only in- 
cluded those objects with the most luminous 
broad lines; objects with weaker broad lines 
may have been missed. Therefore, excluding 
the weak BLRG from the broad line samples 
would have biased BLRG / Q towards higher 
luminosities compared with the NLRG. 

• Weak line radio galaxies: From Figure [5] 
(a) & (b), it is clear that the WLRG lie at the 
low luminosity end of both the MFIR bands. 
Many previous studies have not made the 
distinction between NLRG and WLRG ob- 
jects. In not recognising these radio-loud 
galaxies as a separate population, it is clear 
that the inclusion of WLRG with NLRG 
would lower the mean MFIR luminosity 
of samples of NLRG/WLRG, making the 
BLRG/Q appear relatively str onger emit- 
ters at MFIR wavelengths (see lLaing et~al 



1994 for a a similar discussion relating to 
the [OIII] luminosities). 

Beamed non-thermal MFIR compo- 
nents: When considering the MFIR emis- 
sion from radio-loud galaxies, it is important 
to quantify the degree of non-thermal con- 
tamination of the thermal infrared emission 
from dust. Past investigations that did not 
account for such contamination (e.g. in- 
cluding flat spectrum objects) may have 
suffered a bias in their results. Contam- 
ination from a beamed non-thermal com- 
ponent is more likely in BLRG/Q oriented 
close to the line of sight, and therefore 
could have boosted the MFIR luminosities 
of some of the BLRG/Q objects in pre- 
vious studies. As mentioned in Section 
14.41 we find no more than 24% of our 2Jy 
sample have a possibility of contamination 
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Fig. 9.— Plots of [OIII]A5007 and [OILU3727 vs 5 
GHz radio total radio luminosity. Not including 4 
upper limits in [OIII] and 5 upper limits in [Oil] 

from non-thermal b eamed components (see 



Dicken et all (I2008|)L consisten t with the re- 



sults presented in lClearv et all (|2007l ). 

Although the results of our study support 
the unified schemes, further data are required 
at shorter wavelengths (A < 24//m) to provide 
necessary evidence for greater optical depth in 
NLRG compared with BLRG/Q objects. Addi- 
tionally, our results have been based on the as- 
sumption that the [OIII] emission line luminos- 
ity is isotropic, and it is important to consider 
whether this assumption is valid. 

Previous studies have presented evidence that 
quasars have higher [OIII1A5007 luminosities 
than radio gala xies ( Jackson fc Browne 199(1 
lHaas et al.l l2005lL similar to the results in op- 
position to BLRG/Q and NLRG unification i n 
the MFIR (jHeckman et al.lll994 ISes et al.lll995l) . 
In fact, returning to Figure [21 it is notable that 
there is a weak tendency for the BLRG / Q sources 
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Fig. 10.— Plots of 24^m (top) and 70^m (bottom) 
luminosities normalized by [OIII]A5007 luminosity 
vs 5 GHz radio total radio luminosity. 



to have higher [OIII]A5007 luminosities than the 
NLRG sources fo r a given MFIR l umino sity. On 
the other hand, iHes et alJ (|l993l Il996h demon- 
strate that the [OII]A3727 line, with its lower ion- 
ization potential and critical density, shows no 
evidence for a difference in luminosity between 
BLRG/Q and NLRG. 

In Figure[9]we present L[om]A5007 and £[oii]A3727 
plotted against total radio luminosity (L5GHz). 
A visual inspection of this figure reveals that the 
distributions of the NLRG and BLRG/Q are sim- 
ilar across several orders of magnitudes. Despite 
the BLRG/Q having some of the highest lumi- 
nosities in [OIII] emission lines, there are no sig- 
nificant differences between the two populations 
in thi s plot. This is c o nsistent with the find - 
ings oflLaing et al.l(ll994l) . lMulchaev et al.l(ll994h 
Tadh unter et al.1 ( 19981 ) and Jackson k. Rawlinga 
(|l997h. 



In Figure [10] we again present L 2 4^ and 
plotted against L5GHZ, however in this Figure 
we have normalized the MFIR luminosities by 
L[oni]A5007- If we make the assumption that both 
the [OIII] and MFIR emission are isotropic then 
normalising the MFIR luminosity by £[oin]A5007 
should reveal no difference in the distribution of 
the different optical classifications in the plot. 
However, from a visual inspection of Figure 
riTJl it is clear that, although there is a large 
scatter, the BLRG/Q objects clearly tend to- 
wards lower MFIR/ [OIII] ratios than NLRG; the 
three lowest i24/i/-^[oin]A5007 ratios, and six low- 
est L7o^/i[oin]A5007 ratios are associated with 
BLRG/Q. The BLRG/Q clearly tend to lower 
normalised MFIR luminosities than the NLRG. 
There are three explanations for such a result. 



The [OIII] emission is isotropic, but the 
MFIR emission is weaker for a given [OIII] 
luminosity in the BLRG/Q objects than it 
is in the NLRG objects. This explanation is 
difficult to square with the orientation-based 
unified schemes. In particular, the optical 
[OIII] emission is more likely to suffer dust 
obscuration effects due to the circumnuclear 
torus than the MFIR emission. Moreover, 
assuming that their jets are truly pointing 
closer to the line of sight, the BLRG/Q are 
less likely to be affected by dust obscuration 
than the NLRG. Note, however, that it is 
possible that some individual BLRG / Q may 
have intrinsically under-lumi nous MFIR 
emission (e.g. PKS1932-46: see llnskin et al 
20071) . 



The [OIII] emission is affected by torus dust 
obscuration and therefore appears stronger 
in the BLRG/Q than it is in the NLRG for 
a given isotropic MFIR luminosity. This 
explanation is certainly consistent with the 
orientation-based unified schemes, and the 
conclusions of various other studies (e.g . 
Jackson fc Rawlings 1997, H aas et alj|2005l) . 



It is also consistent with the observational 
evidence for narrow line variability on a 
timescale of years in the B L RG 3C390.3 



(Clavel k Wamstekerl Il987t IZheng et al 



19951 ). which suggests that a proportion of 



the narrow [OIII] emission may be emitted 
on a sufficiently small scale to be affected 
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by torus obscuration. On the basis of Fig- 
ure 10, the degree of anisotropy in [OIII] is a 
factor ~2 - 3 for the majority of the objects, 
only reaching a factor of 10 in a few extreme 
cases. 

• The [OIII] emission is isotropic but the 
Lmfir/ Liqiui ratio decreases with luminos- 
ity because the covering factor of the torus 
decreases as the torus recedes. The reced- 
ing torus model predicts that the BLRG/Q 
objects will have higher [OIII] luminosities 
on average than NLRG, even in the absence 
of dust obscuration, provided that there ex- 
ists a range of intrinsic AGN luminosities 
for a given isotropic radio power. This is 
because the objects hosting the most intrin- 
sically powerful AGN have a greater chance 
of being detected as BLRG/Q. As the inner 
face of the torus recedes with the luminosity 
the covering factor of the torus will decrease. 
This will lead to the following dependence 
of the Lmfir/ L[oiu] ratio on L^i at fixed 

1 /2 

radio power: L M fir/ L[om] °c L b ' ol (See 
Section 16. 3ft assuming a linear relationship 
between £[om] an d i^boi- Therefore in or- 
der to explain the BLRG/Q Lmfir/ -^roni] 
ratios that are a factor w 3x lower than 
the average NLRG, we would require the 
BLRG/Q to be intrinsically more luminous 
than the NLRG at a given radio power by 
a factor of w 10. However as discussed in 
Section 16.31 it is unlikely that this explana- 
tion can apply to the -k70/xm/£[om] ratio 
because the structures emitting the far-IR 
continuum are unlikely to recede as envis- 
aged by the receding torus model. 

7. Summary 

We have presented an analysis of deep Spitzer 
(MIPS) observations for a complete subset of the 
2Jy sample. The main results are as follows: 

1. Heating mechanism: We find tight cor- 
relations between the MFIR luminosities 
and the [OIII]A5007 luminosities, with sim- 
ilar slopes for both 24 and 70/im correla- 
tions (0.75±0.05 at 24^m and 0.72±0.07 
at 70/Ltm). Given that [OIII] is an indica- 
tor of intrinsic AGN strength, we conclude 



that direct AGN illumination is the primary 
heating source for the dust producing both 
the mid-IR and far-IR continuum. These 
correlations are better than those between 
MFIR and radio luminosities, and between 
radio and [OIII] luminosities. 

2. Energetics and dust geometry: Using 
simple arguments we have demonstrated 
that heating of MFIR-emitting dust struc- 
tures by AGN illumination is energetically 
feasible. We identify the dust structure 
producing the mid-IR continuum with the 
circum- nuclear torus, and the dust structure 
producing far-IR continuum with the NLR 
clouds. 

3. Slopes of the correlations: Our finding of 
a non-linear slope in the Lmfir. vs. L^qiii] 
correlation is quantitatively in support of the 
simple receding torus model; however, this 
result can also be explained in terms of a 
non-linear correlation between the strength 
of [OIII] optical emission line luminosity and 
the AGN ionising luminosity. 

4. Unified schemes: Contrary to some pre- 
vious studies we find no strong evidence for 
increased obscuration of MFIR emission of 
NLRG compared to BLRG/Q; the MFIR lu- 
minosities of NLRG and BLRG/Q also cover 
a similar range. These results fully support 
the orientation-based unified schemes, and 
we hypothesise that previous studies may 
have been affected by incomplete, heteroge- 
neous, meager and/or low S/N data. 

5. Weak Line radio galaxies: The WLRG 
in our sample are found to have weak MFIR 
continuum emission as well as weak [OIII] 
emission line emission. This implies that 
the AGN in WLRGs are intrinsically low lu- 
minosity objects, and that their weak op- 
tical emission lines are not simply a con- 
sequence of enhanced dust obscuration by 
circum- nuclear dust. 

6. Starburst heating and triggering: We 

have analysed the increased scatter in the 
distribution of L-jQ^ m vs. £[oni] compared 
to the L24fj,m vs. i[om] correlation. We find 
that the increased scatter is due, in part, to 
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enhanced heating by starbursts in the ob- 
jects in our sample that show evidence for re- 
cent star formation activity at optical wave- 
lengths. The relatively low incidence of ener- 
getically significant starburst activity in our 
sample (15-28%) has implications for our un- 
derstanding of the triggering and evolution 
of radio-loud AGN. In particular, only a mi- 
nority are likely to be triggered close to the 
peak of major gas rich mergers with associ- 
ated starburst activity. Clearly, a variety of 
triggering mechanisms must be present. 

7. Anisotropy of [OIII] emission: Find- 
ing lower ratios of MFIR/[OIII] of BLRG/Q 
compared to those of NLRG, we have con- 
cluded that the most likely explanation for 
this result is that there is mild anisotropy 
(typically factor «2-3) in the [OIII] emis- 
sion due to obscuration by the circumnuclear 
dust in the NLRG that is not seen in the 
BLRG/Q. 

This work is based [in part] on observations 
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is operated by the Jet Propulsion Laboratory, 
California Institute of Technology under a con- 
tract with NASA. This research has made use of 
the NASA/IPAC Extragalactic Database (NED) 
which is operated by the Jet Propulsion Labo- 
ratory, California Institute of Technology, under 
contract with the National Aeronautics and Space 
Administration. Based on observations made with 
ESO Telescopes at the Paranal Observatory. 
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